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We examine the evolution of magnetic fields in an expanding fluid composed of matter and radi- 
ation with particular interest in the evolution of cosmic magnetic fields. We derive the propagation 
velocities and damping rates for relativistic and non-relativistic fast and slow magnetosonic, and 
Alfven waves in the presence of viscous and heat conducting processes. The analysis covers all 
, MHD modes in the radiation diffusion and the free-streaming regimes. When our results are applied 

■ to the evolution of magnetic fields in the early universe, we find that cosmic magnetic fields are 

damped from prior to the epoch of neutrino decoupling up to recombination. Similar to the case of 
sound waves propagating in a demagnetized plasma, fast magnetosonic waves are damped by radi- 
ation diffusion on all scales smaller than the radiation diffusion length. The characteristic damping 
scales are the horizon scale at neutrino decoupling (M u ta 1O _4 M0 in baryons) and the Silk mass 
at recombination (Af 7 « 10 13 Mq in baryons). In contrast, the oscillations of slow magnetosonic 
£ — ' and Alfven waves get overdamped in the radiation diffusion regime, resulting in frozen-in magnetic 

field perturbations. Further damping of these perturbations is possible only if before recombination 
the wave enters a regime in which radiation free-streams on the scale of the perturbation. The 
maximum damping scale of slow magnetosonic and Alfven modes is always smaller than or equal to 
the damping scale of fast magnetosonic waves, and depends on the magnetic field strength and its 
direction relative to the wave vector. 

Our findings have multifold implications for cosmology. The dissipation of magnetic field energy 
, into heat during the epoch of neutrino decoupling ensures that most magnetic field configurations 

generated in the very early universe satisfy big bang nucleosynthesis constraints. Further dissipation 
\Q , before recombination constrains models in which primordial magnetic fields give rise to galactic 

^\ . magnetic fields or density perturbations. Finally, the survival of Alfven and slow magnetosonic 

modes on scales well below the Silk mass may be of significance for the formation of structure on 
'" ^ small scales. 

6 



> 

O 

OO 
O 



I. INTRODUCTION 

In an attempt to explain the origin of galactic magnetic fields through the amplification of primordial fields, several 
authors have considered scenarios for generating magnetic fields in the early universe Q. In such scenarios, one 
attempts to generate fields which will be sufficiently large after recombination at least to seed galactic dynamos and 
at best to produce galactic fields without dynamo amplification. It is generally assumed that after a primordial field 
is generated in the early universe, it becomes frozen into the cosmic plasma and redshifts by flux conservation with 
the expansion of the universe (B cx a~ 2 ; a(t) is the cosmic scale factor). This assumption is usually justified by noting 
that the cosmological plasma is highly conductive and magnetic diffusion is insignificant. 

In this paper, we show that this simple picture of magnetic field evolution is incorrect: at certain epochs in the 
early universe, particularly during recombination and neutrino decoupling, magnetic field energy is converted into 
heat through the damping of magneto-hydrodynamic (MHD) modes. The damping is caused by dissipation in the 
fluid, which arises from the finite mean free path of photons or neutrinos. 

The physical process by which the MHD modes are damped is analogous to that involved in the damping of 
density fluctuations around recombination j^] , and around neutrino decoupling || . Studies of the damping of density 
fluctuations with no magnetic fields present show that, in the diffusive regime (when the scales of interest are much 
larger than the mean free path of photons or neutrinos, ^ m f p ), the effective viscosity and heat conductivity arising from 
the finite mean free path cause the damping of acoustically oscillating density perturbations. Since with the expansion 
of the universe the mean free path of the decoupling particles grows faster than the wavelength of an oscillatory mode, 
all modes whose wavelengths are smaller than the mean free path around decoupling have previously been in the 
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diffusion regime. The rate of damping in this regime ensures that a wave is significantly damped before the mean free 
path of the decoupling particles becomes comparable to the wavelength of the mode. For this reason, the investigation 
of damping in the diffusion regime yields a reasonable estimate of the final damping scales of density fluctuations. 

However, when magnetic fields are added to the fluid, the existence of different MHD modes — Alfven, fast mag- 
netosonic, and slow magnetosonic waves — adds complexity to the problem. We show that while fast magnctosonic 
waves (which include sound waves) damp efficiently in the diffusion regime by the described process, slow magne- 
tosonic and Alfven waves may survive damping by diffusion. Slow magnetosonic and Alfven modes oscillate with 
frequencies which depend on the strength of the background magnetic field and on its direction relative to the mode's 
wave vector, and are in general different from the frequency of sound waves of the same wavelength. In the case of 
a weak background magnetic field or a large angle between the background field and the wave vector, the frequency 
can be slow enough for the damping by viscosity to overcome the oscillation, producing behavior which resembles 
an overdamped oscillator and causes the actual damping of the amplitude to be inefficient. The overdamped slow 
magnetosonic and Alfven modes therefore survive diffusion damping. However, they undergo additional damping if, 
with the expansion of the universe, they enter the so called free-streaming regime, i.e. the mean-free-path grows to be 
much larger that the wavelength of a mode. As a consequence, whereas fast magnetosonic modes are damped mostly 
when radiation is diffusing, slow magnetosonic and Alfven modes are also significantly damped when radiation is 
free-streaming. Therefore, when studying the damping of all MHD modes in order to estimate their damping scales, 
it is necessary to investigate both the free-streaming and the diffusion regimes even before the final stages of the 
decoupling process. 

The damping of MHD modes which causes the dissipation of magnetic energy can be illustrated with the following 
picture: as long as there exist spatially tangled magnetic fields, Lorentz forces accelerate the fluid, setting up oscilla- 
tions about a force-free field configuration; the induced motions are damped by the effective viscosity of photons or 
neutrinos; this causes the exponential decrease in the amplitude of the oscillations and thus results in the straightening 
of magnetic field lines towards a force-free configuration. After the cosmological magnetic fields undergo this damping 
process they have little structure on scales below a characteristic damping scale, and the magnetic energy density in 
such primordial fields is much smaller than that expected from the simple redshift argument above. 

In this paper we follow the evolution of MHD modes and derive their propagation velocities and damping rates 
both in the diffusion and free-streaming regimes during the decoupling of photons and neutrinos. The existence of 
highly relativistic particles with mean free path much shorter than the wavelength of a MHD mode (e.g., photons and 
leptons) requires the use of relativistic MHD. In the radiation diffusion regime, studied in §2, we develop a relativistic 
description of viscous expanding fluids with magnetic fields, while in the free-streaming case, §3, the effects of the 
photons or neutrinos are included through heat exchange and a drag force which they exert on the fluid. This 
procedure allows us to calculate, in §4, the maximum damping lengths after the epochs of neutrino decoupling and 
recombination. Our results may be applied to other astrophysical environments where MHD waves propagate in a 
viscous fluid, since in the derivation of the dispersion relations we leave the sources of viscosity and heat conductivity 
unspecified. 

II. DAMPING OF MAGNETOHYDRODYNAMIC MODES IN THE RADIATION DIFFUSION REGIME 

When the mean free paths of all interacting particle species are shorter than the wavelength of the MHD mode we 
are interested in (A / m fp), it is adequate to study the evolution of a single fluid and account for the effect of the 
diffusing particles by introducing shear viscosity, bulk viscosity, and heat conductivity into the fluid equations |Q| . In 
order to calculate the damping of MHD modes following this approach, we derive linearized relativistic MHD equations 
of an expanding dissipative fluid. We start by reviewing the equations for a non-ideal relativistic fluid in §2.1, and 
add the electromagnetic contributions to the fluid equations in §2.2. In §2.3 we calculate the propagation velocities 
and damping rates for all MHD modes. Our results are applicable for general viscous relativistic and non-relativistic 
plasmas, as long as the pressure is dominated by radiation pressure. 

Throughout the paper we assume that the magnetic field can be decomposed into a large magnitude background 
component Bo(x, t), and a small perturbation, b(x, i). We additionally assume that the curl of the background 
component is negligible when compared to the curl of the perturbations. These two assumptions allow us to solve for 
the damping of MHD modes analytically. 

The use of scalar viscosities and heat conductivity implicitly neglects any anisotropies in these quantities due to 
the presence of the magnetic field. Further, since our equations are derived for an isotropically, homogeneously, and 
adiabatically expanding plasma, the background magnetic field, Bo, is required to have vanishing spatial average on 
sufficiently large scales, (B ) = 0. In our derivation we also neglect gravitational forces because the scales of interest 
are smaller than the Jeans mass scale, and we assume the plasma to be infinitely conducting which is an excellent 
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approximation for most astrophysical plasmas and for the early universe (see, e.g., Ref. 0). 



A. Relativistic Imperfect Fluids 



We consider the evolution of a non-ideal, relativistic fluid in a homogeneously and isotropically expanding back- 



ground using the spatially flat Robertson- Walker metric g^ v = diag(l, 



-a ) and comoving coordinates x^ 



The time dependent scale factor a(t) provides the connection between proper (physical) coordinates x ,tJl and the 
comoving coordinates: x'° = x° and x' 1 = ax z (Greek indices run from to 3 whereas Latin indices run from 1 to 3). 
The relativistic fluid is described by the energy-momentum tensor 



rp\JLU 

EM' 



(1) 



which is separated into three parts: the ideal fluid tensor Tj V , the non-ideal fluid part t^, which accounts for 
dissipation, and the electromagnetic energy-momentum tensor Tg M (added in §2.2). The equations of fluid dynamics 
can be derived from energy-momentum conservation 



0. 



In comoving coordinates Eq. ^| becomes 



dx^ 
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with the dot representing a derivative with respect to time x . 
The energy-momentum tensor for an ideal fluid is 
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where p, p, and U 11 are the total energy density, the total pressure, and the four velocity of the fluid, respectively. 
The non-ideal contributions to the fluid energy-momentum tensor can be written as |q] 
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In this expression, T stands for temperature and "q, £, and k are shear viscosity, bulk viscosity, and heat conductivity 
respectively. 

The effective viscosities and heat conductivity for cither photons or neutrinos are given by 
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(8) 
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where n is the number density of the conserved particles in the fluid and g is the statistical weight of the diffusing 
particles. 
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The exact form of t^ v is partially a matter of definition since, in relativistic fluid mechanics, the fluid velocity 
can be defined either by the flow of conserved particles [|| or by the flow of energy || . These definitions coincide in 
non-relativistic fluid mechanics where the rest mass of particles dominates the total energy. In our case, a relativistic 
one-fluid approximation, the charged and strongly interacting particles (protons, neutrons, electrons, etc.) which 
compose the fluid are all perfectly coupled and have the same velocity as the conserved particle number, the baryon 
number. The energy flow may differ from the particle flow, however, due to the energy transported by the imperfectly 
coupled neutrinos and photons. We choose to follow the velocity of the charged particles (and, therefore, the flow of 
baryon number) which appears explicitly in the magneto-hydrodynamic equations below. 

The conservation of particle number can be written as 







(10) 



where 



nU^ is the particle number four current with n the proper number density of particles. The particle 



number we follow is the net baryon number, n , which is conserved for temperatures below the electroweak transition. 

We can now derive the linearized equations of ordinary relativistic fluid dynamics in an expanding universe from 
Eq. |2| - [l^, by expanding the fluid variables around their background values 



p(x,t) = p (t) 
p(x,t) =p (t) 
T( X ,t) = T (t) 
n b ( X ,t)=n b (t) 



■\-pi(x,t), 
f Pi(x,£), 
+ 7\(x,<), 
+ n$(x,t), 



(11) 
(12) 

(13) 
(14) 
(15) 



The four velocity is that of a stationary fluid element (with respect to the comoving frame) plus a small velocity 
perturbation 



Ctf = (1,0,0,0), u? 



(16) 



We choose in this particular form so that the fluid velocity in proper coordinates, v' = (a/a)x' + v, corresponds 
to an isotropic expansion plus an additional peculiar velocity v. We consider fluids in which the peculiar velocities 
are much smaller than the speed of light, e.g. |v| -C 1, where velocities are measured in units of the speed of light. 
Although the fluid velocities are small, a relativistic treatment is necessary to adequately account for the presence of 
relativistic particles (e.g^, photons and neutrinos). 

Evaluating equations |2| and 10 to lowest order in the fluid variables, we obtain 



£ + »(|W +1 ,)-.e(i 



which represents conservation of entropy when £ = 0, and 



0. 



(17) 



(18) 



represents conservation of baryon number. If we expand the energy momentum tensor and baryon number four-current 
to first order in the perturbation variables pi, p\, T%, and n\, use Eq. and subtract the zeroth order solution, 
we obtain 



dpi 
dt 



Po +Po)-V- v + 3( - ) (p x +pi - 2£V • v) - — V 2 2i 



k d 



(«T V ■ v) = 0, 



(19) 
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and 



3£a 3 dv 



-^V • v + 3 - n? = 0. 
a \a' 



1_ r, _ 9 1^1 



_Vp 1 -4v 2 v--(£ + -7 7 )V(V-v) = 0, (20) 



(21) 



These equations form a complete set describing the evolution of a non-ideal fluid; Eq. [19] represents the first law 
of thermodynamics in local form, Eq. ^0] is the relativistic version of Euler's equation, and Eq. ^l] represents the 
conservation of baryon number. 
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B. Magneto-Hydrodynamics with Dissipation 



We now include the electromagnetic fields. In an inertial frame (denoted by an overhat) the Maxwell tensor has 
the form 



/ E x E y E z 



pia, 



-E x B z 



—E y —B z B x 
V —E z By —B x 



(22) 



where Ei and Bi are the electric and magnetic fields as seen by an observer in the inertial frame. The Maxwell 
tensor in comoving coordinates can be derived from the Maxwell tensor in inertial coordinates (cc M ) by using the 
transformation rules for tensors 



F^ = A» x A\F Xa 



where 



dx^ 



(23) 



(24) 



The coordinate transformation which transforms the locally Minkowski metric — diag(l, — 1, — 1, —1) into the 
Robertson- Walker metric — diag(l, — 1/a 2 , — 1/a 2 , — 1/a 2 ) has 



Thus, in the comoving basis F^ v is 



A"„ = diag(l, 1/a, 1/a, 1/a) 



( E x /a Ey/a E z /a 

-E x /a B z /a 2 -B y /a 2 

-E v /a -B z /a 2 B x /a 2 

\-E z /a By/a 2 -B x /a 2 



The equations of motion for the electromagnetic fields are Maxwell's equations 
and 



dx x dx u dx^ 



d_ 

dx h 



(25) 



(26) 



(27) 



(28) 



where J M is the electric four current. 

In the limit of infinite electrical conductivity, the electric field in the rest frame of the charged particles vanishes 



(29) 
(30) 



This condition, evaluated in the comoving frame using equations 16 and |26j, reads 

E = v x B. 

We decompose the magnetic field into its background value, B , and a small-amplitude perturbation, b(x, t) 

B(x,i) =B (x J *)+b(x,t), (31) 

and impose the following conditions 

b(x,i) <B (x,t) , (32) 



V x B (x,t) «Vx b(x,i) 



(33) 
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to linearize our equations. We can now derive the relevant Maxwell's equations to zeroth and first order in the small 
quantities v and b by using Eq. ^6j ^8], and This yields 



Vb = 0, 

M (fl2b) = ^ V X (V X Bo) ' 
|(a 2 B ) = 0. 



(34) 
(35) 

(36) 



Equation [56] shows that the background field Bo, by flux conservation, redshifts as 1/a 2 with the expansion of the 
fluid. 

To complete the system of equations needed to describe the evolution of the fluid in the presence of electromagnetic 
fields, we must add the contribution from the electromagnetic energy-momentum tensor Tg M to the conservation of 
energy-momentum Eqs. ^ and ^. The energy-momentum tensor for electromagnetism is 



which in comoving coordinates becomes 
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B" 



(37) 



(38) 



(39) 



We can now evaluate the contributions from the electromagnetic stresses to the conservation of entropy Eq. [Tt], 
the first law of thermodynamics in local form Eq. |l9|, and the relativistic version of Euler's equation Eq. ^0|. The 
contribution to the left-hand-side of the zeroth order Eq. [l7] is 



' d ;(a 4 Bl) 



(40) 



87ra 4 dt 

while to first order in v and b, the electromagnetic stresses to be added to the left-hand-sides of Eqs. [l^ and are 
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for Eq. fill and 



1 d 



a 4 b . B ) + (V • v) - - (Bo • V) (v B 



a 4 Bo x (v x Bq) 
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(41) 



(42) 



for Eq. ^0|. After substituting Maxwell's equations (Eq. pi|-p6|) into Eq. fid], we find that the electromagnetic contri- 
bution to Eq. |n] is zero. Similarly Eq. [hi] is identical to zero by virtue of Eq. [56] so that both Eq. [l?] and |l9] are 
unmodified. The only coupling between the field and the fluid to first order occurs through the velocity of charged 
particles and the curl of the magnetic field. 

Note that the first term of Eq. fi^2] is only important in the relativistic limit. For a mode with frequency lu and 
wavenumber k, Eq. [3^ implies tub ~ kvB. Therefore, the relative contribution of the first term in Eq. ^ compared to 
the second term is of order (oj/k) 2 . Hence, the first term can only be neglected when the group velocity of a mode, 
duj/dk sa w/k, is much smaller than the speed of light. 

To first order in the quantities Ti, n\, pi, Pi, v, b, equations [l9|, ^l], [5~4"|-[36"|, together with the equation obtained by 
adding Eq. ^ to the left-hand-side of Eq. describe magneto-hydrodynamics in an expanding fluid. The following 
definitions help to rewrite our equations into a more convenient form 



6=£ 
T 



(43) 
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r, B ^ ?» r , (44) 

(47r(p r +p r )) 2 (4:n{p r +p r )) 2 

V = V , £' = ~ . «' = KT , (45) 

3(pr+Pr) ' 3(p r +p r ) ' ' 3(p r +Pr) 

where p r — gi^T^ /30 and p r = p r /3 are the average energy density and pressure of relativistic particles, and g is the 
total statistical weight of relativistic particles. 

For a fluid comprised of baryons and relativistic particles (e.g., photons, neutrinos, e ± -pairs, etc.), the energy 
density and pressure up to first order in the small quantities are given by 

P = Po + Pi 

= Pr (l + 45)+p b (l + A), (46) 
V = Pa + Pi 

= p r (l + AS), (47) 



where p^ = nJmN is the baryon energy density, and is the nucleon rest mass. In writing Eq. 47, we assume that 
baryonic pressure is negligible in comparison to radiation pressure. In this case, and when £ = 0, Eq. [l^ and Eq. [l8| 
imply simple redshift relations for the temperature To ~ 1/a and the baryon number density n\ ~ 1/a 3 . We also 
define 

R(t) ee (48) 

Ap r 

as a measure of the relative importance of baryon mass density with respect to energy density in relativistic particles. 
For R — > 0, both energy density and pressure are dominated by relativistic particles, whereas for R 3> 1 the energy 
density is dominated by the baryon rest mass and the pressure is dominated by radiation. 
In terms of the newly defined variables, the equations of magneto-hydrodynamics become 

a 1 w n ( ®' \ £' n K 'w2r k! d . 

a 



<5 + 36£'- S+— V-v-6 - )-V - v - —V 2 5 — V-v = 0, (49) 

\a J 3a \a J a ar a at 

2x 



I + ^ (!) ) fa + R > - ^ + >) - ( v2v + l v < v • v 



^V(V-v)+B x ^xB j+^Box(Vxb)=0, (50) 



^b= -Vx(vx B ) , (51) 
at a v ' 



V • b = , (52) 



|B = . (53) 



C. Dispersion Relations 

In order to calculate propagation velocities and damping rates, we first derive the dispersion relations for the 
different MHD modes by Fourier transforming all perturbative variables (5, A, b, and v, generically represented by 
$ below) using the convention 

$(x,i) =y"d 3 k$(k,t)exp(ik-x) (54) 
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in which k is a constant comoving wave vector. 

The time dependence of $(k, t) is modified by the expansion of the fluid which introduces a time variation into the 
frequency and the amplitude of the modes. With this in mind it is convenient to write 



$(k, t) = $ k (t) exp [ I iw(t) dtj . (55) 

The decrease of the amplitude due to damping is included in the exponential part through imaginary solutions for uj 
while the explicit time dependence of $k(*) acounts only for the effects of the expansion. 

The system of equations resulting from the substitution of Eq. into Eq. 49 ~ 53| is solved separately for the different 



MHD modes: Alfven waves, for which the density and the temperature of the fluid are uniform and the velocity of 
the fluid is perpendicular to the background magnetic field (k = fcx,B = B x x + B z z,b = by, v = vy and 5 = 0); 
and magnetosonic waves, for which the velocity of the fluid makes an arbitrary angle with the background field 
(k = fcx, B = B x x + B y y, b = by, v = t) x x + v y y). Note that sound waves, which propagate along the background 
field without affecting it (Bo || k | v and b = 0), are a special case of magnetosonic waves. In all dispersion relations 
and their solutions, 9 denotes the angle between the background magnetic field and the wave vector. 

The dispersion relation for magnetosonic modes allows two solutions: slow magnetosonic modes and fast magne- 
tosonic modes. Fast magnetosonic modes are similar in nature to sound waves, while slow magnetosonic modes are 
closer in character to Alfven waves. This fact plays an important role in the damping of magnetosonic waves, and 
is apparent for weak magnetic fields (pb <C pfluid) where fast magnetosonic modes oscillate with v almost along the 
direction of k and involve oscillating density perturbations, while slow magnetosonic modes oscillate almost perpen- 
dicularly to k and have close to vanishing density perturbations. In the special case k || B fast magnetosonic waves 
become sound waves and there are no slow magnetosonic solutions. (For a discussion of MHD modes see, for example, 
Ref. §. 



We obtain the dispersion relations by substituting Eq. 55 into Eqs. 49H52L The dispersion relations are derived 



to first order in k', rf ' , which corresponds to the lowest non-trivial expansion in powers of ^ m f p /A (where A is the 
wavelength of a mode). We use the WKB approximation neglecting the time derivatives of the Fourier amplitudes, 
d$k(t)/ dt <C ui, and the time derivative of the frequency, dui/dt -C —uj 2 , which arises from the d 2 v/dt 2 term in 
Eq. |50|. This approximation is valid for modes with oscillation frequencies much higher than the expansion rate, 
uj > H. 

This procedure for Alfven waves yields: 

3uj 3 k' + iuj 2 (l + R + B^ + uj(^R+3r)'(^J ^ -iB 2 cos 2 e(^j =0. (56) 

Complete dispersion relations for all MHD modes are, due to their length, placed in Appendix [X| 

In the following two sections, we present the solutions for oscillation frequencies and damping rates derived from 
the dispersion relations in the following two limits: the oscillatory limit, when the solution is oscillatory with an 
exponentially decaying amplitude; and the overdamped limit, when the amplitude of modes exponentially decrease 
without completing an oscillation. 



1. Oscillatory Limit 

The solutions to the dispersion relations, uj, generally consist of a real and an imaginary part, which represent the 
oscillation frequency and the damping rate, respectively. In the oscillatory limit, the dissipative effects are such that 
the fluid oscillates many times as it damps, Rew <C Imu. In this case, the dispersion relations can be solved by 
considering all the viscosity and heat conductivity terms as perturbations on the ideal fluid dispersion relation. 

The solutions to the dispersion relations will be conveniently expressed in terms of the speed of sound: 



dp\ 1 



dpj s y/3(l + R) ' 



(57) 



and the relativistic Alfven speed 

VA — — 

'l + R + B 2 



VA = ; B ° ■ (58) 
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The relativistic Alfven speed includes the magnetic field energy density in the denominator, which ensures that for 
strong magnetic fields the Alfven speed does not exceed the speed of light. 

For all modes in the oscillatory limit we first solve the dispersion relations for an ideal fluid and then compute the 
first order contributions from the dissipative terms. For clarity, all the solutions in their general form have been placed 
in Appendix [A|. In this and the following section we give the solutions for each MHD mode in the cosmologically 
relevant limit of weak magnetic fields (Bq -C 1) and negligible redshift terms (a/a). 

For weak magnetic fields, the leading terms in the frequencies for fast magnetosonic waves do not depend either 
on the magnetic field strength or on the direction of propagation, and are therefore the same as the frequencies for 
sound waves: 

"™ = + 'MV + TTr^ + WTR) ( ') (:)' + ot6l) - (59) 

This reproduces the solution for propagation and damping of sound waves given in Weinberg (1971). Similarly, the 
frequencies of slow magnetosonic waves and Alfven waves are identical to leading order in Bq and have the following 
form: 

The solutions show that, while for small magnetic fields the damping of slow magnetosonic and Alfven waves 
proceeds through shear viscosity, fast magnetosonic waves are damped by shear and bulk viscosity, as well as heat 
conductivity. Furthermore, fast magnetosonic waves damp differently in different regimes: they damp predominantly 
by heat conductivity when the matter density is larger than the radiation density, and by shear viscosity when the 
radiation density dominates. Note that after taking the time dependence of all variables into account, the expansion 
of the fluid affects the frequencies directly through a/ a terms (see Appendix |A|), and indirectly through the integral 
in Eq. |55|. For instance, as in the case of sound waves when R <C 1 , the oscillation frequency of a fast magnetosonic 
mode with a given wavelength in a radiation dominated expansion is twice the frequency of the same wavelength 
mode in a static background metric. 



2. Overdamped Limit 

When dissipative effects become very strong, oscillations of MHD modes are inhibited and the evolution of a given 
MHD mode is dominated by the exponential decay of its amplitude with time. We seek solutions in the extremely 
overdamped regime by expanding the equations in powers of Re w sc/Imaj sc, where lo osc is the frequency of a wave 
derived in the oscillatory limit. 

In general, a dispersion relation expanded in powers of Rew osc /Imw osc has several solutions distinguished in nature 
by their initial conditions. For example, in the case of Alfven waves, fast decaying solutions arise from initial 
conditions such that when the velocities of the fluid are damped away by shear viscosities the amplitude of the 
magnetic perturbations vanish as well. In contrast, when initial conditions generate slowly decaying Alfven modes, 
the fluctuations are not erased as the velocities damp to zero; after the damping of fluid motions, the remaining 
magnetic forces tend to accelerate the fluid, although inefficiently because of the strong viscous damping. Since 
energy dissipation rates are proportional to the peculiar fluid velocity, the timescale for dissipation of the magnetic 
field perturbation of slow decaying modes may be extremely large. 

While the amplitudes of the fast decaying modes damp at rates similar to the ones calculated in the oscillatory 
regime (§2.3.1), the amplitudes of slowly decaying modes decay at significantly different rates. For weak magnetic 
fields (Bq <C 1), the decay rate for the amplitude of overdamped slow magnetosonic modes is 



,SM _ • /„,2 / 'M , ■ A 



2 „,2 „„2 



V a COS 



and the decay rate for Alfven modes is 



%k'v%[-\ + z ^ , (61) 

^od = « A 3r] , ■ (62) 

Note that all modes with relativistic propagation velocities (Rew osc ~ k/a) never enter the overdamped regime in 
the diffusion limit. For this reason a discussion of overdamped relativistic sound and fast magnetosonic waves is not 
necessary. 
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III. DAMPING OF MAGNETOHYDRODYNAMIC MODES IN THE RADIATION FREE-STREAMING 

LIMIT 



Slow magnetosonic and Alfven modes which become overdamped during the diffusion regime survive the damping 
and with the expansion of the universe enter the free-streaming regime when the mean free path of the decoupling 
particles grows to be larger than the wavelength of a mode. In order to investigate the additional damping that 
these modes undergo in the free-streaming regime, we study the general case of MHD in an expanding fluid in the 
presence of a uniform background. Similar to our analysis of the diffusion regime, we study the evolution of a single 
dissipative fluid. However, in this case the fluid is comprised of all the particles with mean free paths much shorter 
than the wavelength of the MHD mode, while the decoupling particle species, whose mean free path is now large as it 
decouples from the rest of the fluid, represents a uniform background on the scales of interest. The dissipation arises 
from occasional collisions of the fluid particles with the relativistic background. 

We generically define a drag coefficient a and a heat exchange coefficient 7 in the following way: the drag force per 
unit volume on the fluid element from scattering with the background particles is given by 



(63) 



and the heat exchanged between the fluid element and the background is 

C?Pthcrmal _ T\ 

77. = iTrT Pthcrmal • 

ot T 



(64) 



The exact form of these coefficients is obtained by calculating the transfer of momentum and heat per scattering and 
averaging it over the distribution of background and fluid particles. We presented the coefficients later separately for 
neutrino decoupling and photon decoupling. 

In order to derive the free-streaming fluid equations, we use the fluid equations developed in §2, as well as techniques 
for finding solutions described therein. The heating rate from Eq. 64 is incorporated into Eq. njl while the drag 
force from Eq. ^ is added to Eq. ^ together with the magnetic field contribution from Eq. ]42[ Although the 
local thermodynamic equilibrium between the fluid and the free-streaming component does not hold in general, we 
consider the case in which the mean scattering time between particles of the fluid component and the free-streaming 
component is shorter than the characteristic expansion time scale of the fluid. In this case the temperature and 
velocity of the background is the same as the average temperature and velocity of the fluid. All other assumptions, 
including non-relativistic fluid velocities, are carried over from §2. 

The resulting equations are: 



dpi 
dt 



1 



(po +Po)V 



v + 3-(pi 
a 



Pi) 



f m Pt hernial 
^0 



(65) 



i__a 

^dt 



1 1 dv 

{a 4 (Po +Po)v) + -Vf>i + — B x (-J 
a 4-7T v ot 



x B n 



i-Box(Vxb) 



-avpo 



(66) 



ot a a 







(67) 



51 



Eq. 53 they 



and together with Maxwell's equations Eq. 

To derive dispersion relations for a given fluid, we have to 
and substitute these into the above set of equations. This 
baryonic fluid with free-streaming photons; and a fluid which 
and e + e~ pairs, in a background of free-streaming neutrinos 
fast magnetosonic modes are given in the appendices. Here 
magnetosonic modes and Alfven modes, the two modes that 
free-streaming regime before recombination. 



form a complete set. 

specify the energy density, matter density, and pressure, 
is done in the rest of §3 for two fluid combinations: a 
consists of baryons and relativistic particles like photons 
All the dispersion relations as well as the solutions for 
we present the solutions to dispersion relations for slow 
in the presence of weak magnetic fields survive into the 



A. Neutrino Free-Streaming Limit 



Around neutrino decoupling the fluid consists of tightly coupled baryons, photons, and e e pairs. The dominant 
component of the pressure is radiation pressure 
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and, since R ~ 0, the speed of sound is 

- = 75 • (69) 

The energy density (equal to pthormai) has contributions from all relativistic particles, counted in the number of degrees 
of freedom g r , 

Pr=9r^T i . (70) 

The heat exchange coefficient and drag coefficient defined by Eq. ^34] and Eq. |6^ may be computed by averaging the 
transfer of energy in each scattering between a distribution of background particles and a distribution of fluid particles 
JiT| . For neutrino free-streaming they have the form: 

7 ~ a w n w — = — — , (71) 

and 

a ~ 7 , (72) 

where the second equation follows on dinensional grounds. Here a w is the cross section for scattering of neutrinos 
with other weakly interacting particles, n w is the number density of weakly interacting particles (scatterers), g v is the 
neutrino statistical weigth, and l v is the neutrino mean free path. 

Following the steps used in the diffusion regime (§2.3), we obtain dispersion relations for the different MHD modes 
and present them in Appendix pL 



1. Oscillatory Limit 

The oscillation frequencies and damping rates for slow magnetosonic and Alfven modes are again obtained by first 
solving the dispersion relations for an ideal fluid and then solving for the first order dissipative terms. In terms of 
the previously defined Alfven speed, the solutions for small magnetic fields are the same for slow magnetosonic and 
Alfven modes and have the form: 

^ A = ±^cos^)+fa. (73) 
The frequencies for oscillatory fast magnetosonic waves (including sound waves) are presented in Appendix H. 



2. Overdamped Limit 

As in the diffusion regime, the solutions in the extremely overdamped regime (Reu osc <C Imu; osc ) are derived by 
expanding the equations in powers of Rew osc /Imw osc . The overdamped solution in the case of weak magnetic fields 
is, for slow magnetosonic modes 

SM i 2 Aiv\ cos 2 ( k\ 2 

and for Alfven waves: 

A Aiv\ cos 2 6 ( k 



3 a 



(75) 



As in the case of radiation diffusion (§2.3), relativistic fast magnetosonic modes do not become overdamped in the 
free-streaming regime. 
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B. Photon Free-Streaming Limit 



Fluid equations for the modes in the free-streaming limit around photon decoupling are somewhat different from 
the cases analyzed so far. Namely, the only contribution to the energy density of the perturbations is the thermal 
energy density of baryons, 

3 

^thermal = ~ («e + U p )T = 3n fc T, (76) 

which enters Eq. |65| but can be neglected in Eq. ^36] because it is much smaller than the matter density. Here n e and n p 
are electron and proton number densities respectively. Furthermore, since the photons can be considered decoupled 
on the free-streaming scales, the only pressure left to support the oscillations is the pressure of the baryonic fluid 
itself: 

Pb = (n e + n p )T = 2n b T. (77) 

This yields the gradient of pressure in Eq. ^ which depends both on density and temperature fluctuations and is best 
expressed through the sound speed: 

ivp = ~c s 2 (VA + V<J), (78) 
P 5 

where c s is the adiabatic baryonic speed of sound for a fully ionized proton-electron fluid 



opb ) s y 3 m p 

With these substitutions for the densities and the pressure, we obtain dispersion relations for different MHD modes 
and present them in Appendix 

The drag and heat exchange coefficients, which appear in the dispersion relations and their solutions, are similarly 
obtained as in the neutrino free-streaming case, and have the following form jyj : 

a ~ a T n e — = • r —, (80) 

Pb '7 -ft 



and 



7~— ^a. (81) 

m e 



1. Oscillatory Limit 

Unlike in the relativistic cases when the photon pressure dominates, the structure of the non-relativistic equations 
with free-streaming photons allows for oscillating magnetosonic modes with two different propagation velocities and 
damping rates. These modes are commonly referred to as adiabatic and isothermal, depending if heat transport is 
slow or rapid compared to the oscillation time: a mode is adiabatic when to 3> 7 and u > a, and it is isothermal 
when 7 > w > a. Alfven modes have only one solution since they do not include density or temperature fluctuations 
and therefore are not affected by heat transport. 

Again we derive the dispersion relations from Eq. |5l|-[53|and Eq. [35|-|67|and place it in Appendix^. Before presenting 
the solutions for slow magnetosonic and Alfven modes, it is useful to introduce the non-relativistic Alfven speed, 



The oscillation frequency and the damping rate for slow magnetosonic waves in the adiabatic regime TLeu ^> 7, a 
are given by: 

cu „ f k\ (a 7 B n sin 2 9 \ 
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and 



= ±C S COS 



where the upper solution corresponds to weak magnetic fields such that c s 3> ca and the lower solution to strong 
magnetic fields with c s <c^. It is important to remember that both of these solutions are derived for a background 
magnetic field whose energy density is much smaller than the energy density in photons. The condition for adiabaticity 
is dependent on the strength of the magnetic field since oscillation frequencies of magnetosonic waves are different 
for strong and weak magnetic fields. Slow magnetosonic modes in the isothermal regime, 7 > Reu > a, have the 
solution 

sm ■ „fk\ .a 

u 

and 



u™ = ±c A caa0[ -)+%-, (85) 



^ = ±f-c s cose(^j+q, (86) 

where again the upper solution is for c s 3> ca and the lower solution is for c s <C ca- 
The result for Alfven waves is 

u>i c = ±c A cosd(J\+i^. (87) 

The frequency of non-relativistic fast magnetosonic waves in the adiabatic and isothermal limits are placed in 
Appendix O. 



2. Overdamped Limit 

The frequencies for slow magnetosonic waves corresponding to slow exponential decay are: 



sm _ • c a cos2 9 ( k 



2 

a , 



and 

SM . 3cj cos 2 9 ( k 

u od = 1 ; — 

5a 



(89) 



with the upper solution for ca <C c s , and the lower solution for ca ^> c s . The slowly decaying mode of overdamped 
Alfven waves has frequency 



a c 4 cos 2 9 ( k 



(90) 



In contrast to relativistic MHD, some non-relativistic fast magnetosonic modes enter the overdamped regime. Their 
damping rates can also be found in Appendix |c[ 



IV. DAMPING OF MAGNETIC FIELDS IN THE EARLY UNIVERSE 



In this section we discuss the implications of the damping of magneto-hydrodynamic modes for the evolution of 
cosmological magnetic fields. Magnetic fields generated in the early universe are likely to be randomly oriented, 
spatially varying fields with small coherence lengths, usually of the order of the horizon at the epoch when the fields 
were created. We assume that the magnetic fields are created with magnetic field energy below equipartition with the 
radiation energy density, i.e. Bo ^ 1- For an arbitrary magnetic field configuration, we choose a separation of scales 
such that in a given volume the field can be described as an approximately force-free background magnetic field Bo, 
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and a spectrum of propagating modes b(k), where |b| <C |Bo | - In this case we can decompose the propagating modes 
into slow and fast magnetosonic, and Alfven modes with different wave vectors k and different phases. Although the 
condition |b| <C |Bo| may not be easily achieved for every field configuration, the predicted evolution of propagating 
modes is indicative of the general field evolution. In particular, the efficient viscous damping discussed in this paper 
should cause the dissipation of magnetic energy in generic field configurations. 

We are interested in the evolution of individual MHD modes from before the epoch of neutrino decoupling to 
recombination. For each epoch we wish to determine the characteristic scales over which pre-existing cosmic magnetic 
fields are damped. As previously described, the evolution of fast magnetosonic waves is distinctively different from the 
evolution of slow magnetosonic modes and Alfven waves; therefore they are discussed separately: fast magnetosonic 
waves in §4.1, and slow magnetosonic and Alfven waves in §4.2. 

For the calculation of the damping scales we need the expressions for the mean free path of the decoupling particles 
as well as the ratio of the baryon density to the photon density. While the universe cools from temperatures below 
the electroweak breaking scale (T ~ 100 GeV) to neutrino decoupling (T„ ~ MeV), neutrinos are the particles with 
the longest mean free path and therefore the most efficient momentum and heat transporters. The neutrino mean 
free path at temperature T can be written as 

where ni and n q are lepton and quark number densities, Gf is Fermi's constant, and G F T 2 is a typical weak interaction 
cross section. The quantities gi and g q are the statistical weights of relativistic weakly interacting leptons and quarks 
present at the epoch of interest. 

At temperatures below the completion of the e + e~-annihilation (T ~ 20 keV), heat and momentum are most 
efficiently transported by photons. The dominant process limiting the photon mean free path during this period is 
Thomson scattering of photons off electrons which gives the photon mean free path the following form: 

w * j_ * 10 -c m (^y _> (^)"W' . w> 

where X e is the number of free electrons per baryon, gt is the Thomson cross section, and n e is the electron density. 

The damping of MHD waves is particularly efficient during the epochs of neutrino decoupling and recombination, 
when efficient momentum transfer and heat transport arise from the quickly growing mean free path of the decoupling 
particles. Therefore, all the variables in the above equations have been scaled to their values at those epochs. Also, 
since at neutrino decoupling the dominant scattering process is scattering of neutrinos off leptons, the appropriate 
values for the parameters in Eq. |9l] are gi = 8.75 (e^ and six neutrino species) and g q = 0. The remaining unspecified 
parameter, the ionization fraction X e , drops within a short time from 1 to ~ 10~ 5 during recombination which occurs 
approximately at Ti? ~ 0.25 e V. 

The baryon mass density is negligible when compared to the radiation energy density during neutrino decoupling 
(R ~ 0), while around recombination it is approximately given by 

^^^VY^Y (93) 



4p r \0.25eVy V0-0125/ 

where is the fractional contribution of baryons to the closure density and h is the present Hubble constant in units 



of 100 km s 1 Mpc . In writing Eq. 93 we implicitly assume that neutrinos have decoupled from the remaining 
particle species. 



A. Damping of Fast Magnetosonic Waves in the Early Universe 



The damping of all fast magnetosonic modes is to leading order equivalent to the damping of sound waves if the 
energy density in the large-scale magnetic field is much smaller than the energy density in radiation. The damping 
occurs due to the diffusion of either neutrinos prior to neutrino decoupling, or photons before recombination. From the 
expressions in Eqs. ^ and |9^, which represent mean free paths in proper units, it can be seen that in both cases the 
comoving mean free path grows with the expansion of the universe. As a consequence, MHD modes with wavelengths 
X p = 2Tra(T)/k smaller than the mean free path l m fp{T) at time t(T) were in the diffusion regime at some prior time 
in the early universe, i.e. A P (T") > Z m f p (T') at t(T') . 
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The amount of damping that fast magnetosonic modes undergo in the diffusion regime can be calculated using the 
damping rates in Eq. [59]. From the leading damping term, 2ir]' (k / a) 2 / (1 + R), and the definitions in Eqs. and [7], 
we see that the amplitude of the perturbation is damped between time t = and t by 



exp - ^dt> . (94) 



t J 

fc mfp 

Af ' 

If we define a characteristic damping scale as the largest comoving wavelength of an MHD mode whose initial amplitude 
has been damped by at least one e-folding by time t, this damping scale approximately corresponds to the comoving 
photon/neutrino diffusion length, 

,'if J '' < 95 > 

which is the distance a photon/neutrino has random walked between time t = and t. 

In this section we present all the damping scales in comoving units, in particular, comoving to the present epoch, un- 
less stated otherwise. The diffusion damping scale calculated for fast magnetosonic waves prior to neutrino decoupling 
(T > 1 MeV) is: 

A™ ~ 2 x 10 20 cm ( -*-) " f ( \ l V 9 -l±°lY 1 . (96) 

\MeV J \10.75J \5.25j \ 8.75 J V ; 

Note that the damping scale at neutrino decoupling converted to proper units approximately corresponds to the causal 
horizon at that time, A„(MeV) « 5 x 10 10 cm. 

In a similar fashion we can compute the comoving damping scale of fast magnetosonic waves due to the effects of 
the finite photon mean free path at lower temperatures. Around recombination, where we assume matter domination, 
this damping scale is 

A™~7x 10 25 cmf— ^— ^ 4 (^L) 2 (Sl Q h 2 Y k . (97) 
7 V°- 25eV / V 0.0125/ V ' V ' 

In this expression £l is the total density in units of the critical density at the present epoch. 
Baryonic mass scales can be associated with the damping length scales by defining 

where pb is the average baryon mass density at temperature T. The baryonic mass scale associated with the damping 
scale of fast magnetosonic waves by neutrino diffusion around neutrino decoupling is 

M r - 1O" 4 M (^-) ^ (^-Y 1 (J*-) f ( ?L±9lY 1 (ShOL) , ( 99 ) 
°VMeVy \10.75J \5.25J \ 8.75 J V.0-0125 / ' V ' 

and by photon diffusion around recombination 



^ = ^tolfrj (SJHS) < 100) 

The above calculations for the damping scales are approximate in that we assume that the damped modes are in 
the diffusion regime, l vn <ti X p . The diffusion approximation is not a valid approximation late within the decoupling 
epochs. Around neutrino decoupling, while our analysis uncovers the order of magnitude of the damping mass scale, 
an improvement on the diffusion approximation only, would not yield more accurate results since the calculated 
damping scale approximately corresponds to the causal horizon at that time. Around recombination, our treatment is 
analogous to the calculations of the damping of sound waves in Ref . [Q ; in writing Eq. |9?] we assume an instantaneous 
recombination while the ionization fraction X e decreases gradually to zero during recombination. More detailed 
treatments for the damping of sound waves during recombination have been performed using the Boltzman equation 
[p~2|— |i~4|| or using a two-fluid model f p^yi^ ]. Since the dominant damping terms of fast magnetosonic waves are the 
same as those of sound waves, values for the sound wave damping scale calculated in more accurate models may be 
used as better approximations to the fast magnetosonic damping scale. A review of previous Silk scale calculations is 
given in Ref. 
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B. Damping of Slow Magnetosonic and Alfven Waves in the Early Universe 



Unlike the damping of fast magnetosonic waves, the damping of slow magnetosonic and Alfven waves in the early 
universe proceeds through several different stages. We illustrate these stages by following an Alfven wave from after 
neutrino decoupling up to recombination. 

Initially, in the diffusion regime where A 3> Z m fp, a wave oscillates and damps in the same manner as described 
for fast magnetosonic waves. The oscillation frequency and the damping rate are shown in Eq. and after using 
Eqs. ^ and ^ they approximately become 

k\ „ . /fc x J 



w^~±VA\-JcoaO + ilj^-J . (101) 

The damping rate in this expression is the same as the damping rate for fast magnetosonic modes, and is valid as 
long as vaCOsO 3> l^k/a. The crucial difference, however, is that for a small background magnetic field in the early 
universe, the oscillation frequency of an Alfven mode (u^fc/a) is much smaller than the oscillation frequency of a 
fast magnetosonic mode with the same wavelength (v s k/a). While all fast magnetosonic modes of interest satisfy 
the condition for damping in the oscillatory regime (v s 3> Lfe/a), an Alfven mode can become overdamped when, 
with the expansion of the universe, the mean-free-path becomes large enough for dissipative effects to overcome the 
oscillation (vacosQ ~ l 7 k/a). One may define a temperature dependent wavelength A oc j, 

X od (T,B )^^l, (102) 
va cos 6 

such that modes with proper wavelength A p > A 0( j are oscillatory while modes with A p < A ot j are overdamped. 

Overdamped modes are a superposition of fast and slowly decaying overdamped modes, and the relative amplitudes 
depend on the phase of the mode when it becomes overdamped. Fast decaying overdamped modes damp at rates 
similar to the oscillatory modes, and therefore their damping is equivalent to the damping of fast magnetosonic modes 
discussed in the previous section. In this section, we follow the significantly different evolution of the slowly decaying 
overdamped modes which experience the least damping. The amplitude of a slowly decaying overdamped Alfven 
mode damps with a rate given by Eq. ^2], which in terms of the photon mean-free-path is 

A .v 2 A cos 2 

- • (103) 

Since the damping rate is inversely proportional to the growing mean free path, the integrated damping rate is much 
smaller than the integrated damping rate of fast magnetosonic modes during the same period. As a result, the 
damping in overdamped diffusion is inefficient and the damping scales of Alfven modes at the end of the diffusion 
regime are smaller than the damping scales of fast magnetosonic modes. 

The left-hand-side of Fig. [I] (left of the dotted line) illustrates the evolution of Alfven waves in the diffusion regime 
for cos6> = 1 and a background magnetic field of Bo w 10~ 3 (corresponding to Bq w 3 x 10~ 9 Gauss at present). In 
the diffusion regime there are two important temperature dependent scales: the photon diffusion length scale (dashed 
line), which is the scale damped by one e- folding by the time the universe cools to temperature T provided that modes 
are still in the oscillatory regime; and the overdamping length scale given by Eq. |102| (solid line), which shows the 
temperature at which a mode with comoving wavelength A c enters the overdamped regime. The modes which do not 
damp significantly in the radiation diffusion regime before they become overdamped, preserve their amplitude until 
they reach the free-streaming regime, apart from small additional damping during the transition itself. On the graph, 
these are all the modes with comoving wavelength larger than that given by the intersection of the solid line with the 
dashed line. Therefore, the intersection of these two lines roughly represents the largest Alfven mode that is damped 
by one e-folding in the diffusion regime. Its position depends on the strength of the background magnetic field and 
on the angle between the field and the wave vector. 

Some overdamped modes enter the free-streaming regime before recombination (to the right of the dotted line in 
Fig. |l|) as the mean free path of the decoupling particles grows with the expansion. The dissipation coefficients in 
the free-streaming regime are inversely proportional to the mean free path unlike those in the diffusion regime. This 
implies that, when a wave enters the free-streaming regime, it is initially overdamped and becomes oscillatory when 
the drag force and the heat conduction decrease as the mean free path increases. It also implies that modes in the 
free-streaming regime undergo most damping while overdamped. The damping rates during free-streaming, derived 
using Eq. BO in Eqs. and B(X are: 
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when a wave is overdamped, and 



when it oscillates. 



•( g) (£) «>s^. (104) 



w^sicAf-^cosff + ti^L. (105) 



Using Eq. 104 in Eq. B3, we find that for Alfven waves during overdamped free-streaming, the largest comoving 



wavelength whose amplitude is damped by one e-folding at time t is: 



(A A (/)r [\ 2 A cos 2 9 R(t>) !^dt' = c 2 A cos 2 9 d 2 (t)R(t) , (106) 



where d 7 (i) is the comoving photon diffusion length. As we can see, the damping depends on the strength of the 
background magnetic field and the angle between the field and the wave vector. The above damping length evaluated 
at recombination is the characteristic Alfven wave damping length for the free-streaming regime, because all modes 
stay overdamped before recombination regardless of the strength of the magnetic field. 

The damping of Alfven waves with free-streaming photons is illustrated on the right-hand side of Fig. [l]. The 



dashed line shows the damping scale from Eq. 106 . The solid line marks t he tr ansition from overdamped to oscillatory 
behavior, defined in the same fashion as in the diffusion regime (Eq. |102| ), with the overdamped region to the 
left of the line. Since all the modes that cross the dashed line before recombination get damped by one e-folding 
during overdamped free-streaming, the length scale marked by this line at recombination represents the free-streaming 
damping scale for Alfven modes. 

In our analysis we have assumed the WKB approximation. This approximation does not formally hold during all 
the discussed epochs for every mode. In particular, the approximation breaks down for comoving wavelengths between 
the damping length scale for Alfven waves in the photon free-streaming regime, A 7 , and the photon mean free path 
(i.e., between the dashed and dotted lines on the right-hand-side of Fig. ||). Although the damping rates calculated 
in the preceding section predict no damping in this region, some damping is in principle possible since here WKB is 
not a good approximation. In Fig. |^ we present the result of a numerical integration of the magnetic field amplitude 
of Alfven waves in different wavelengths through the epoch of interest. The mode with the largest wavelength which 
still damps before the epoch of recombination (T ~ 0.25eV) is consistent with our analytic estimate. Note that in 
this calculation we assume that recombination never occurs (i.e., X e ~ 1 for T < 0.25eV in Fig. ||). We have extended 
the calculation into this non-physical regime to learn when the results obtained by the WKB approximation deviate 
significantly from the numerical results. From the plot we see that significant discrepancies happen only with modes 
that would damp after recombination if X e was kept fixed. 

The damping of slow magnetosonic waves proceeds similarly to the illustrated damping of Alfven waves except for 
two differences. First, slow magnetosonic waves in overdamped diffusion damp at a slightly different rate than Alfven 
waves because of the additional damping expressed through the extra term in Eq. ^l|. This damping is relevant only 
for the waves that propagate at very large angles with respect to the background magnetic field. Second, the damping 
rates for slow magnetosonic waves during free-streaming depend on whether the non-relativistic sound speed is larger 
or smaller than the non-relativistic Alfven velocity (Eq. [38|-Eq. |S9| ) . The damping scale when c s > ca is the same as 
the damping scale for Alfven waves, while the damping scale for c s < ca is unique to slow magnetosonic modes. 

We can define damping length and mass scales for Alfven and slow magnetosonic waves analogously to the previ- 
ous section, and consider that the waves below these scales would have dissipated by the time of neutrino/photon 
decoupling. The general evolution of all the damping length scales with temperature, followed through the different 
damping stages up to decoupling, can be found in Appendix [d| In this section we present the final damping lengths 
at neutrino decoupling (T ~ lMeV) and at recombination (T ~ 0.25cV), for the background magnetic field below 
3 x 10~ 9 Gauss, which is the current estimate of the upper limit on magnetic fields on Mpc scales (see, e.g. Ref. [p7[). 

The final damping scale of Alfven waves for this range of the background magnetic field strengths is determined 
by the free-streaming damping length. At neutrino decoupling, most modes in free-streaming are damped while 
overdamped, although the largest modes are damped just as they begin to oscillate again. This determines the 
comoving damping scale at neutrino decoupling to be: 

M'~B.„.(^y(&) i (e±±) i . dor, 
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The background magnetic field is expressed through Bg = Bo/ (3 x 10~ 9 Gauss), where Bo is the background magnetic 
field strength scaled to the present epoch. During photon decoupling, the largest Alfven mode damped by one e-folding 
during free-streaming is damped while still in the overdamped regime, and its comoving wavelength is: 



A^~2x lO 23 cm-B 9 cos0( Q ti 



n b h 2 

0.0125 



(108) 



The baryonic mass scales which correspond to the damping length scales at neutrino and photon decoupling are 

,2 n / s -3 



M A ^ 1O" 13 M B% cos 



0.0125 / V 10.75 



9r 



9v_ 
5.25 



9i + gg 
8.75 



(109) 



and 



M// ~ 10 6 M Q 5|cos 3 0| <>„//" 



n b h 2 v 

0.0125 ) 



(110) 



The characteristic damping length sca le for slow magnetosonic waves at neutrino decoupling is the same as the 
damping length of Alfven waves (Eq. 107 ). At recombination however, the slow magnetosonic damping scale depends 
on the value of the non-relativistic sound speed relative to the non-relativistic Alfven speed. When c s > ca the 
damping scale for slow magnetosonic waves in the free-streaming regime is the same as the damping scale for Alfven 
waves (Eq. |l08j). On the other hand, when c s < ca the slow magnetosonic damping scale is: 



A™ ~ 3 x 10 21 cm cos ( n h' 



(111) 



which is independent of the background magnetic field strength. The condition c s > ca is equivalent to Bq ^ 
5 x 10 -11 (r2fc/i 2 /0.0125) 1 / 2 Gauss for the strength of the large scale magnetic field at the present epoch, or Bq < 
2 x 10~ 5 (n b h 2 /0.0125)^ 2 and B 9 < 1.7 x 10~ 2 (fi 6 /i 2 /0.0125) 1 / 2 for the two different scalings of the magnetic field 
strength present in our equations. 

The baryonic mass scales associated with the damping lengths of slow magnetosonic modes are: 



Mr. 



10 



" 13 M 0J B|cos 3 i 



and 



0.0125 / V 10.75 



fjr 
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5.25 



9i +9g 
8.75 



(112) 



M. 
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1O 6 M Bg cos 3 9[ tt h 2 



1M Q cos 3 9 ( n h 2 



0.0125 / 



n b h 2 \ 

0.0125 J 



for c s > ca 
for c s < c A 



(113) 



V. CONCLUSIONS 



In this paper we have studied the effects of dissipation on the propagation of MHD modes in an expanding fluid 
composed of matter and radiation. We have derived the propagation velocities and damping rates for fast and slow 
magnetosonic, and Alfven waves in the radiation diffusion and radiation free-streaming regimes. The derived damping 
rates have general applications in magnetized relativistic and nonrelativistic astrophysical environments. We have 
applied the damping rates to the evolution of MHD modes in the early universe to show that cosmic magnetic fields 
suffer significant damping from before neutrino decoupling to the end of recombination. 

Fast magnetosonic waves are damped by radiation diffusion on all scales smaller than the radiation diffusion length 
in analogy to the propagation of sound waves in a demagnetized plasma. The characteristic damping scales are: 
the horizon scale at neutrino decoupling, M™ s» 10~ 4 Af Q in baryons, and the Silk mass at recombination, M^ M « 
10 13 M Q in baryons. In contrast to fast magnetosonic waves, slow magnetosonic and Alfven waves reach an overdamped 
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regime during which the damping is not very efficient; further significant damping occurs once the radiation is free- 
streaming on the scale of the perturbation. The maximum damping scales for slow magnetosonic and Alfven modes 
in general depend on the magnetic field strength and the direction of propagation with respect to the background 
magnetic field. At neutrino decoupling the damping scale is M^- SM 10~ 13 MqBq cos 3 9, the same for both types of 
modes. At recombination, if B 9 < 1.7 x 1CT 2 , this scale is M^' SM m 10 6 M Q B| cos 3 9, and if B 9 > 1.7 x 1(T 2 , the 
damping scale for slow magnetosonic modes is different: M^ M « 1M@ cos 3 9. The background magnetic field strength 
scaled to the present epoch is expressed in terms of B 9 — (Bo/3 x 10~ 9 Gauss), the current observational limit on the 
large scale magnetic field. 

Our findings have various implications for cosmological magnetic fields. The dissipation of magnetic field energy 
into heat during the epoch of neutrino decoupling ensures that most magnetic field configurations generated prior to 
neutrino decoupling satisfy big bang nucleosynthesis constraints. The observed element abundances require that the 
energy density in magnetic fields be less than one-third of the photon energy density during nucleosynthesis |p^| , |l9| . 
Even if processes prior to neutrino decoupling generate magnetic fields with initial energy density comparable to the 
photon energy density, neutrino damping causes the magnetic energy density to decrease substantially relative to that 
of radiation by the time of nucleosynthesis. 

Further dissipation before recombination lessens considerably the magnetic field energy in primordial magnetic 
fields available for generating galactic magnetic fields or density perturbations. The damping of fast magnetosonic 
modes implies that a sizable fraction of the magnetic energy density is erased up to the Silk scale. 

Finally, although Alfven and slow magnetosonic modes also undergo significant damping, their damping scales 
depend on the strength and the direction of the background magnetic field and are generally smaller than the damping 
scale for fast magnetosonic modes. As long as mode coupling is not effective, which we expect to be true for non- 
oscillating modes before recombination, magnetic energy density can be stored in Alfven and slow modes on scales 
well below the Silk mass. The survival of these modes may be of significance to the formation of structure on relatively 
small scales. In particular, these modes may be responsible for fragmentation of early structures as well as seeding 
early star formation. We will further discuss the consequences of the cosmological evolution of MHD modes in a 
subsequent paper. 
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FIG. 1. Scales relevant for the evolution of Alfven and slow magnetosonic waves before recombination, calculated for modes 
that propagate at cos# = 1 and a background magnetic field of Bo = 3 x 10 -9 Gauss today. All length scales are given in 
comoving units. Any mode with fixed comoving wavelength A c will at cosmic temperature T be in the photon diffusion regime 
if it is to the left of the dotted line, or in the photon free-streaming regime if it is to the right. Modes with wavelength A c will 
at temperature T be non-oscillatory (overdamped) if they are between the two solid lines. The two dashed lines indicate the 
temperature at which a mode of given wavelength is damped by one e-fold, either during its oscillatory evolution in photon 
diffusion or its overdamped evolution in photon free-streaming. The figure assumes Qt = 0.0125 and h = 1, and equality 
between radiation and matter energy density at Teq = 5.5eV. See §4.2 for a more detailed explanation. 



FIG. 2. The evolution of the Fourier amplitude as a function of cosmic temperature calculated numerically (solid line), and 
analytically using the WKB approximation (dashed line), for Alfven waves with three different comoving wavelengths indicated 
on the figure. For the calculation we assume Bo = 3 x 10 _9 Gauss, cost9 = 1, fif, = 0.0125, and h = 1, and we have fixed the 
ionization fraction at X e — 1 even for temperatures below the approximate temperature of recombination T ~ 0.25eV. The 
largest scale which the numerical calculation shows to be damped by one e-fold before recombination (at ), corresponds to the 
analytically calculated length scale of 2 x 10 23 cm. 



APPENDIX A: RADIATION DIFFUSION 



This appendix contains the diffusion regime dispersion relations and their solutions. The dispersion relation for 
magnetosonic waves, expanded to first order in k' , ij' , and £' (which corresponds to a first order expansion in l 7 /X), 
is: 

lu 5 (6k'(1 + R) + 3k'B%) 

,4/7i , D\2 , 52/ 



2 / 7 \ 2 / 7 \ 2 



+ ita A ({l + Rf + B 2 (l + i?)) 

! (77/(1 + R) (?) + 3£'(1 + R) (?) ' + k'{R 2 

+ V 'B 2 cos 2 o(?) 2 + 3£'B 2 cos 2 of?) - 2k'B 2 cos 2 $(?) * + 2^R(1 + R) + -^R)) 
-3C'SoCos 2 0[ - ) + 2k' Bl cos 2 1 - ] - Z B'qR 



3 a \ a J a 

+ s^i^f +M i^ +M ^. l , (A1) 

Since the scales of interest are smaller than the horizon scale, we only keep a/a terms which are first order in l^/du 
and A/dn- In addition, note that there are two terms in the last row which seem to be second order in the expansion 
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variables, although for some slow magnetosonic waves for which J3ocos# is so small that Bq cos 2 6 <C BqT]'k' or 
Bq cos 2 6 -c B^a'a/a, these terms play a crucial role in determining the damping rates for overdamped solutions. 

We recover the dispersion relation for sound waves from the one for magnetosonic waves, if Bq = 0, or if the 
magnetosonic wave propagates along the field lines, k|| Bq : 



3wV + iu> z (l + R)+uj 



The dispersion relation for Alfven waves is: 



-5 ; = ° 



3wV + iuj 2 I 1 + R + Bl ) + oj I -i? + 377' f - 



sii^cos 2 6>( - ) =0 



(A2) 



(A3) 



The solutions to the dispersion relations are presented using the following convention: the oscillatory part of the 
frequency is denoted by coq = Rew, and the damping rate by u>\ = Imu. For Alfven waves the solutions are: 



LUq = + 



1 



(1 + R + B 2 ) 
1 

2(l + R + B 2 ) 



^B 2 (l + R + B 2 )[-\ --\~\ i? 2 , 



1 /d 
4V a 



3r/ - +3k'^cos 2 6» - 



(l + R + Bi) 



a -R 

a 



(A4) 



For magnetosonic waves, a general solution for a non-dissipative fluid is: 



v s I k 



UJQ = ± 



1 + g + 3£ 2 (1 + g) + 5 2 cos 2 g 
{l + R + Bi) 
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[l + R + W$(l + i?) + i?o cos2 
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12£^cos 2 6>- 



1 + 



1 + -R + Bn) 



(A5) 



This solution contains two magnetosonic solutions: fast magnetosonic, whose fre que ncy is obtained by taking the plus 
sign, and slow magnetosonic, using the minus sign. For weak magnetic fields Eq. A5 may be expanded to second order 
in Bq, and the solutions become: 
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The imaginary (dissipative) parts of the magnetosonic frequencies, again to first order in a/a, k', and rj', and to 
second order in Bq are: 



, FM 



and 
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1 2(1 + R)\i 
For sound waves Eq. |A8| reduces to: 
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which reproduces the solution of Weinberg (1971) if the expansion is neglected (a/a = and a = 1). 



(A8) 



(A9) 



(A10) 



(All) 
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APPENDIX B: NEUTRINO FREE-STREAMING 



During neutrino free-streaming, the dispersion relations are: 

uj 2 (1 + Bl) - -iau - Bl cos 2 6»( - I =0 



(Bl) 



for Alfven waves, and 



lo^I + BD-iu 3 - 
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Bn — 
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Bl cos 2 6 fk 
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= 



(B2) 



for magnetosonic waves. 

The solution to the dispersion relation for oscillatory Alfven waves is: 



, v\ cos 6 ( k\ 3 a 
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The oscillatory magnetosonic solutions are expanded to first order in the dissipation coefficients a and 7, corre- 
sponding to the first order expansion in \jl v . This yields for fast magnetosonic waves: 
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and for slow magnetosonic waves: 



,SM 



±u A cos 9(1 - Bl sin 2 9) ( - ] + \i(a + 3B 2 7 ) 
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The decay rates for overdamped slow manetosonic and Alfven waves are given in the text in Eq. |7| and Eq. ||. 



APPENDIX C: PHOTON FREE-STREAMING 



During photon free-streaming, the dispersion relation for Alfven waves is: 



to — no I a H — 
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Bl cos 2 6» fk 
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and the dispersion relation for magnetosonic waves is 

'k 



w°-iu*[ 2a + 7 + 3- l ■„/' 



c 2 | - 



2 n 2 . fir 



~RVa' +" 2 + 2a 7 + 2-(2a + 7 ) 



,2 . B 2 \fk 



R \a 



^ 2 1 ^0 N / & 



R I \ a 



a 



a 2 I 7 



20:7- 



'fc 
■ cj 1 — 



i? U/ 'V5 s i? / al 5 s R I a R 



| i7 ^(i,.„. (C2) 



The magnetosonic dispersion relation has been derived to first order in a/a, but without any approximation in the 
dissipative coefficients a and 7. 

The oscillatory solution to the dispersion relation for Alfven waves is: 



22 



osc 



k\ if a 
a I 2 \ a 



(C3) 



For fast magnetosonic waves in the adiabatic regime (Rew 3> 7, a) the oscillatory solution is: 



(C4) 



2B% J\aJ \2 5 2 a 

and in the isothermal regime (7 S> Rew 3> a) it is: 
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The condition for adiabaticity is dependent on the strength of the magnetic field when c A c s . 
In the limit Bq = 0, the above solutions give the solutions for sound waves: 

«=W*W?+;+*V (06) 

in the adiabatic regime (c s k/a 3> 7, a) and 



a / V 2 5 a 



in the isothermal regime (7»c s t/a>a). 

The oscillation frequency and the damping rate for slow magnetosonic waves in the adiabatic regime are: 
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and in the isothermal regime: 
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Note that all the solutions are derived for a background magnetic field whose energy density is much smaller than the 
energy density in photons (Bq <C 1). 

In contrast to relativistic MHD there do exist non-relativistic, fast magnetosonic waves in the overdamped limit. 
The decay rates for the amplitudes of these overdamped fast magnetosonic waves are: 

(CIO) 

.fir for 

The decay rates of overdamped slow magnetosonic and Alfven waves are given in the text in Eq. - Eq. pu. 
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APPENDIX D: DAMPING SCALES 



In this appendix we give the evolution of the damping scale as a function of temperature for Alfven and slow 
magnetosonic modes in the early universe. The temperature dependence of the damping scale of fast magnctosonic 
waves is given in Section 4.1. Before neutrino decoupling, the damping scale for Alfven and slow magnetosonic modes 
evolves as follows: 
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The damping scale at T > T\ is approximately the same as the damping scale for fast magnetosonic modes, and 
represents the diffusion length of the decoupling particles. It is illustrated for photon decoupling by the dashed line 
on the left-hand-side of Fig. |l|. The largest wavelength mode still damped by one e-fold during oscillatory diffusion 
becomes overdamped at temperature T\\ 
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i.e. the temperature at which the oscillatory damping scale (diffusion length) and the overdamping length on the 
graph intersect. 

Since the damping is inefficient during the overdamped diffusion phase, the largest length scale damped in diffusion 
represents the maximum damping scale, until further damping during free-streaming. Temperature T 2 represents the 
point when the overdamped free-streaming damping scale (equivalent to the dashed line on the right-hand-side of 
Fig. [I]) exceeds the maximum scale damped so far, 



= io (iwr 
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If this happens before decoupling, the damping scale grows further, and is now determined by the free-streaming 
damping scale. Finally, the largest length scale damped in free-streaming is the one which still damps by one e-fold 
before becoming oscillatory again, determined by the intersection of the dashed and the solid line on the right-hand- 
side of Fig. pi). This happens at: 
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The final damping scale is determined by this transition; as it is given in Eq. Dl for T < T3, it includes the scales 
damped during the transition itself. Note that for large magnetic field strength, with energy density approaching 
equipartition with radiation energy density (Bg ss 10 3 ), and cos# « 1, the maximum damping length scale is deter- 
mined by damping during the oscillatory diffusion phase rather than during the overdamped free-streaming phase. 
In this case the damping scale of Alfven and slow magnetosonic modes becomes similar to the damping scale of fast 
magnetosonic modes. 

Around photon decoupling, the damping scale for Alfven modes, and the damping scale for slow magnetosonic 
modes when c s > ca, have the same form: 
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On the other hand, if ca > c s , the damping scale for slow magnetosonic modes evolves like: 
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The transition temperatures are determined by matter-radiation equality at Teq = 5.5 eV(i7o/i 2 ), and by the last 
scale damped by one e-fold in oscillatory diffusion, which becomes ovcrdamped at Ti, 
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and presents the largest damping scale until damping in free-streaming damps even larger scales, at T < T 2 with 
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However, if T 2 < Teq the damping scale of 10 21 cm (B 9 cos 6») 3 (^ 6 /i 2 /0.0125) is valid until (T/0.25eV) 



-1/4 



2TJ g 7 5 cos6'- 8 / 5 (r2 6 /i 2 /0.0125)~ 4 / 5 (r2n^ 2 /0.0125)- 1 / 5 when it changes to 3 x 10 21 cm cos 0(T/O.25eV)- 5 / 4 (fi o ^ 2 ) -1 

The above equations are derived for the case in which the universe is radiation dominated during damping in 
oscillatory diffusion (Ti > Teq). This is always true for the background magnetic fields with B$ < 1. Although for 
larger background magnetic fields the damping processes are the same, the temperature dependence of the diffusion 
damping scale and its time of overdamping might be different, resulting in a different final damping scale. For example, 
for BgcosO w 10 3 (or vacos# ~ 1), the oscillatory diffusion damping scale dominates up to recombination, and its 
evolution in the matter dominated era determines the damping scale at recombination. 

Slow magnetosonic modes which propagate almost perpendicular to the background magnetic field have different 
damping scales due to the additional term in Eq. pTl However, this scale is substantially different from those given 

8 / i • 

in the text only for modes with cos# < 0.03_B g , which makes the influence of this additional damping term on the 
overall damping of slow magnetosonic modes negligible. 
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